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Abstract 



By using a usual instanton method we obtain the energy splitting due to 
quantum tunneling through the triple well barrier. It is shown that the term 
related to the midpoint of the energy splitting in propagator is quite different 
from that of double well case, in that it is proportional to the algebraic average 
of the frequencies of the left and central wells. 
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I. INTRODUCTION 



The Euclidean time formalism of path integral has been applied successfully to the prob- 
lems in which quantum tunneling occurs [|l],|2| . In most cases, their results are consistent with 
the WKB approximation || as well as numerical calculations Q. Recently much attention 
has been paid to periodic instanton method which can be applied to high-energy quan- 
tum tunneling and has a possibility to explore the baryon- and lepton-number violation in 
electroweak theory. Furthermore, there is another approach called valley instanton method 
that allows exploration of the quantum tunneling from false to true vacuua 

In general the Euclidean time formalism gives the information about the ground state, 
i.e., the ground state wavefunction and the ground state energy. In the case of double well, 
it allows applications of path integral by opening a possiblity of finding a classical path in 
Euclidean space, and explains quantum tunneling through the barrier between wells result- 
ing in an energy splitting of the degenerate ground states. This energy splitting was also 
investigated through the WKB approximation and numerical calculation as well as through 
the path integral with Euclidean time (instanton method). These results are consistent with 
each other and it was conjectured that the instanton method may give more exact result 
than the WKB approximation [3]. 

We here consider a triple well potential on the ground of the instanton method. For 
simplicity, we take the triple well potential in the following form; 



This potential has three wells among which two side wells have topologically identical vacuua, 
and central well is rather broader as shown in Fig. 1. This potential form may be realized 
by designing supperlattices with different compound compositions whose applications might 
come true in electronic devices such as laser without inversion || . Therefore, understanding 
the structure of the energy level, which influences the dynamics of electrons, is essential in 
constructing electronic devices. 

Compared with the double well case, there are two different aspects. Firstly, the energy 
spectrum forms discrete blocks, each of which is composed of three levels. The level located 
in the middle of each block has different symmetry; for example, energy levels in the lowest 
block is arranged in the order of even, odd and even states in energy. Secondly, the central 
vacuum is not identical with the others locating at the sides. This indicates that the structure 
of the block can be adjusted by changing the curvature of the central vacuum. This aspect 
might be important in practical electronic engineering. 

In this paper we obtain the central position and width of the lowest block, i.e., energy 
splitting, and partial information about the wavefunctions of the lowest block using the 
instanton method. The term related to the midpoint of energy splitting turns out to be 
proportional to algebraic average of the frequencies of the left and central wells. In Sec. II, 
we find out what physical quantities are related to the propagator with large Euclidean time. 
One instanton contribution to the propagator is calculated in Sec. Ill, and in order to obtain 
a desired form of the propagator we use the dilute gas approximation for the configurations 
composed of multi-instantons and anti-instantons in Sec. IV. In Sec.V a brief conclusion is 
given. Throughout this paper, we take % = 1 and the mass m — 1 for simplicity. 
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II. PROPAGATOR AT LARGE EUCLIDEAN TIME 



We consider a triple well problem in which barriers are so high that at least one discrete 
eigenvalue can exist below the barrier height. Generally in the case of triple well potential, 
an energy level is split into three levels with slightly different energy values due to the 
quantum tunneling through barriers. 

Let \L >, | Co >, and \Ro > denote the normalized lowest energy eigenstates of the 
isolated left, central, and right well, repectively. When barriers are appropriately high and 
wide, the three lowest eigenstates in triple well can be obtained as following through a 



variational method 10 



\E' > ~ . 1 (|L > +a + \C > +\R >) (2.1) 



2 + a 2 + 

\E[> ~-L(|L >-|i? >) (2.2) 

\E' 2 > ~ , 1 (\L > -a_|C > +\Rq >) (2.3) 
a?_ 



where 



1 




(2.4) 



Here, 



H LL =< L \H\L >, H cc =< C Q \H\C >, H LC =< L \H\C >, (2.5) 

where H is the total Hamiltonian. In getting these states, we use the fact that the ground 
state wavefunctions of each well are taken to be almost separated due to large potential 
barriers and there is a two-fold symmetry. The lowest three eigenstates in Eqs. (|2.1| - |2.3|) are 
orthogonal to each other and have following eigenvalues: 



< E' \H\E' >= E' ~ ^LL^u + iicc^ ^+\^lc\ (2 g) 



2H LL + a 2 + H C c ~ 4:a + \H LC \ 
2 + a2 

<Si|ff|^>=^~ff ii (2.7) 

< E>\H\E> E> * ^ LL + alHcc + ^\H LC \ _ ^ 

If the curvatures of the potential with three vacuua are equivalent, then this problem is 
reduced to three-fold degenerate one. In this case, H LL = H C c and a± = y/2. Thus the 
lowest energy eigenstate is split into three eigenstates having equal spacing AE = \Z2\Hlc\ 
between the eigenvalues. 

It is well known that for a large Euclidean time interval the propagator gives the infor- 
mation of the ground state energy. In the triple well case, the propagator from left(</>j) to 
central vacuum(0/) is approximately described as follows: 

< 0/|e- 2OT |^ >~£< <P f \E' 3 >< E'^ > e- 2E 'i T . (2.9) 

j=0 
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Here, since the wave function of the state \E[ > is odd , < 4>f\E[ > vanishes, and the 
coefficients of the terms corresponding to j = and j = 2 in the Eq. (|2.9| ) are equal, that is, 



2a+ 2a_ 



2 + a\ 2 + a 2 + s 

Then, the propagator becomes 

-2HT I j, 2 xiJn U ^ „-2ST, 



(2.10) 



< 0/|e- 2iii |0, >~ < ^|C >< L o |0 4 > e" 2 ^ sinh2A£T, (2.11) 



where 



E= 2 , AE = ^ ^. (2.12) 



This result shows that, in the case of triple well, the propagator for large Euclidean time T 
tells us about the width and the central position of the lowest energy block. 

In the following sections we will calculate the propagator by using the instanton method, 
and show that the propagator with the dilute instanton gas approximation takes the same 
form as Eq. ( p.ll|) . 



III. ONE-INSTANTON CONTRIBUTION 

In the present section we will calculate the one-instanton contribution to the propagator. 
The Euclidean propagator from (0 i; — T) to (0/,T) is written as 

<<P f \e~ 2HT \<p i >= [** T D<f>e- S °M, (3.1) 

J d>i,-T 



where the Euclidean action is 



S E [cj>] = i T dr 

J-T 



\<P 2 + Vtf) 



(3.2) 



The path giving major contribution is determined by the classical equation of motion, i.e., 

5Se/5<J) = 0. For the triple well potential in Eq.(l.l) whose central well is wider than the 
side wells, the equation of motion becomes 

= = 2a0(30 4 - 4/? 2 2 + /? 4 ). (3.3) 

The classical solution of this equation of motion, so-called instanton, is 

»d(r) = i == ■ (3-4) 

<l _|_ e 2/32Vto(r- to) 



We here take a boundary condition as (pi = —f3 and <f)f = (from left to central vacuum) 
Substituting the solution into Eq. ( p.2|) , and taking a T — > oo limit, we obtain 
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lim S E [M = ^ • (3.5) 

In order to evaluate the contribution of paths neighboring on the instanton solution to 
the propagator, we define 

<P(r) = <f )d (r)+v(r), (3.6) 

where t](t) is an arbitrary function satisfying a boundary condition: 

V (-T) = V (T) = 0. (3.7) 

Now the action can be divided into the two parts: one from instanton and the other from 
the paths neighboring on the instanton. Therefore, the total action becomes 



T 

S F , = S E [4> d } + I dr 

-T 



hf + aY{r) V 2 



(3.8) 



where 



Y(t) = 1504(r) - 12/3 2 2 (r) + p\ (3.9) 
In Eq. (|3.8|) , we keep up to the second order terms of 77. Then the propagator can be written 



as 



< 0\e~ 2HT \ -(3>= e- SE[ ^ ] I(T), (3.10) 



where 



and 



I(T) = / Drie' J dTvMri , (3.11) 

J(-T,0) 



M=-\^- 2 +aY(r). (3.12) 

Our next task is to calculate the functional integral I(T) which has approximated to a 
Gaussian form. This is, however, not simple because it contains a zero mode, which is easily 
seen from the fact, 

M0d(r) = O. (3.13) 

It indicates that <p c i (r) is the eigenvector corresponding to zero eigenvalue of M. Therefore, 
the zero mode should be treated carefully in calculating I(T). Here we follow the method 
used by Gildener and Patrascioiu || in performing the functional integral. 

Let the eigenfunctions and corresponding eigenvalues of M be ipm and e m , respectively, 

Ml\) m = Emlpm- (3.14) 

Then the small fluctuating path rj (r) can be expanded as a linear combination of the eigen- 
functions: 
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V( T ) = ^C m 1p m (T). 



(3.15) 



In the Gaussian functional integral I(r), the integral variable rj can now be changed to the 
coefficients c m , which gives a formal solution including the Jacobian in the form 



dc n 



dc ' 5 \l e ' 

UL -n m= Q V tm 



(3.16) 



To surmount the difficulty of the zero mode, one can use the collective coordinate method, 
which gives 



where 



I(T) = 2T^[0 d ]/ o , 

° c n m= i V £ m 



(3.17) 



(3.18) 



In order to find an explicit expression for I we use the method of changing variables 

gr 

I(T) = \ 2itN(T)N(-T) 



ll| . Then the functional integral I{T) can be represented as 

T dr 



where 



N{t) = <p cl 



1/2 



— T N 2 (t) 



'I _|_ e 2f3 2 V2ary3/2 



(3.19) 



(3.20) 



For large T, this becomes 



/(T) = /3 (^v^VV^V^t 



Comparing Eq. (|3.21| ) with Eq. (|3.16|) and Eq. (|3.18|) , we obtain 



Jo = ^(2/5 2 v / 2^£ ) 1/2 e^ /32v/ ^ T . 



(3.21) 



(3.22) 



Imposing the boundary condition, ip n (—T) = ip n (T) = 0, yields the lowest energy e in 
WKB approximation at large T[3] as 



s = 8aP 4 e~ 2 ^ T . 
Inserting this into Eq. ( |3.22| ), we get 



4/5 3 



7T 



(3.23) 



(3.24) 
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From Eq. Q3.17D , therefore, the final form of the functional integral I(T) is expressed in the 
form 

I{T) = -P 3 TjsM(aV2^) 1/2 e-l^ T . (3.25) 

7T v - 

In the case of double well, it is well known that the exponent in the exponential function 
is closely related to the curvature of vacuum, that is, if the second derivative of the double 
well potential at vacuum positions is u 2 , then the exponential function is given by e~ wT . 
However, in our triple well potential the curvature of the side vacuua is different from that 
of the central vacuum (Fig.l). Note that in this case the exponent is also related to the 
potential curvatures at the vacuum positions in a way of algebraic averaging, that is, 

- - = \^ (3.26) 

where 

uj\ = V"{-(5) = 8a/? 4 , uj\ = V"(0) = 2a/? 4 . (3.27) 

As will be seen shortly, u/2 is the central position of the energy block composed of three 
almost degenerate lowest levels. 

From Eq.( |3.10| ), the contribution of one instanton to the propagator becomes 



< 0|e- 2 " T | -/3>= 2KT^S E [M ] J^- s ^e^ T , (3.28) 

where 

n = A^J^. (3.29) 
V on 

In the next section the contribution of multi-instanton to the propagator will be explored. 



IV. MULTI-INSTANTON CONTRIBUTION 

In the previous section we consider the contribution of the stationary solution(instanton) 
and its neighboring paths. However, since the path integral in Eq. ( |3.1| ) includes all possible 
paths starting from (—(3, —T) and ending at (0,T), let us consider the paths consisting of 
rii instantons and n<i anti-instantons which are located at large distance with each other in 
Euclidean time so that the interaction between them can be ignored. In order to satisfy the 
boundary condition, the number of instanton must be one larger than that of anti-instanton. 
Furthermore, since the total number of possible configurations is 2™ 2 , the multi-instanton 
contribution to the propagator is 

<0|e- 2HT |-/3> 

oo oo 9™2 , rrj 

= E E o ^-n 2 ,i (ni + n2) , e-(" 1+ "^^(2KT v /^[0d]) rai+n y- l 
,U.,1J , 

to (2n + 1 



e wT 



wJsmJ^.-™ t p«- M "yyW , (4 , ; 
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Using the identity, 



00 x n 1 ^ 

S(^^ smhVi ' (4 - 2) 



the infinite sum in Eq.(fOj) is represented by a hyperbolic sine function. Therefore, the final 
result is of the following form, 



< 0|e~ 2HT | -P>= J^-e-^ smh(2V2KTJS E [(j) d ]e- SE[ ^ ] ). (4.3) 

V Z7T V 

Comparing this result with Eq. fl2.1l|) , we get the central position E and the width 2AE of 
the lowest energy block in terms of potential parameters, a and (5, which are 

E = ^[3 2 V2^ (4.4) 

and 

AE = i/^(2«) 3 / 4 /5 4 e-^ 4 ^. (4.5) 

V 37T 

For sufficiently large (5, in other words, both large well separation and high potential barrier, 
the energy splitting decreases exponentially as expected. Moreover, the product of the 
amplitudes of the ground state wavefunction at vacuum positions is expressed as 

P(3 —x 1 / 2 



<0\E' ><E' \-P>=^-V2a) • (4.6) 



V. CONCLUSION 

The quantum tunneling for the triple well potential case is investigated by means of the 
instanton method. Unlike the case of double well potential, the curvature of the central well 
is topologically different from those of the others, which results in the fact that the midpoint 
of energy splitting is proportional to the algebraic average of the frequencies of the left and 
central wells. We also obtained the energy splitting AE of the lowest block. Although this 
kind of approach can be extended to the n-tuple(n > 3) well, much complication in counting 
the possible configuration gets involved. In the case of quadruple well, for example, two 
different types of instanton and anti-instanton solutions should be considered carefully. The 
study about this subject is in progress. 

ACKNOWLED GEMENT 
This work was supported in part by Nondirected Research Fund, Korea Research 
Foundation, 1995 and by Korea Science and Engineering Foundation (961-0201-005-1). 



8 



REFERENCES 



[1] S.Coleman, "The uses of Instantons" , in The why's of subnuclear physics, ed. A.Zichichi 
(Plenum, New York, 1979). 

R.Rajaraman, Solitons and Instantons, (North- Holland, the Netherlands, 1987). 
E.Gildener and A.Patrascioiu, Phys. Rev. D16, 423 (1977). 
K.Banerjee and S.P.Bhatnagar, Phys. Rev. D18, 4767 (1978). 
S.Yu.Khlebnikov, V.A.Rubakov, and P.G.Tinyakov, Nucl. Phys. B367, 334 (1991). 
J.Q.Liang and H.J.W.Muller-Kirsten, Phys. Rev. D50, 6519 (1994); 51, 718 (1995). 
H.Aoyama, T.Harano, M.Sato, and S.Wada, Nucl. Phys. B466, 127 (1996). 
H.Aoyama, T.Harano, H.Kikuchi, M.Sato, and S.Wada, [hep-th/9606159| . 
S.E.Harris, Phys. Rev. Lett. 62, 1033 (1989). 

R.L.Liboff, Introductory Quantum Mechanics, (Addison- Wesley Publishing Company, 
Inc. 1992). 

[11] R.F.Dashen, B.Hasslacher, and A.Neveu, Phys. Rev. D10, 4114 (1974). 



9 



FIGURES 

FIG. 1. Triple well potential as a function of position. 
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